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Pee OO VERVIEW 

The effectiveness of 2 system depends not only on its 
ability to meet its specified perfornance requirements, but 
meron on its ability to perform when needed, for the duration 
of its assigned missions, and for its operational lifetins. 
The technical disciplines concerned with «hese time-related 
Seem Characteristics are reliability, maintainabilicy and 
ieogescic Support. 

Commonly used methods for prediction and démonstration 
of corrective mainténance repair times assume the validity 
Seecne lognormal distribution. 

Mamcrae> to help in f2cuSing aztention on these matters, 
a statistical analysis on data sets of demcns*rated and 

i2ld repair times has be2n conducted as part of a study on 
miewapplication of the lognormal and gamma distributions <o 


q 


COTrrective maintenance data. The results of the analysis 
Meemgi.ven in this thesis. A previous study of alectronic 


Systems and equipments concluded shat the lognormal dis-ri- 


te 


rh 


Mieeeon Wal a suitabl= descriptor for repair times [Re 


Bee FS URPOSE AND APPROACH 


The objectives sf this study have been; 





Goo fo test the 


(wD 
(D 


xponential, lognormal, and gamma 
Meer ibutions as descriptors of repair time data on mechani- 
cal and electronic systems and equipments, 

fp to Vendeverua: one Jdeoga2='ual or gamma distri - 
bution is a suitable descriptor for repair times, 

(c) to estimat= the percentage error caused in 
assuming €xponentiality for availability and maintainability 


Meeeulat>oOns when in fact the distribution is lognormal. 


2. Systems and Data Analyzed 


Perihny-Six Sets Of repair cime data for slectronic 


and mechanical systems / 2qQuipments wer2 analyz 
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Date £Or the electronic systen (Sets 8-20) cam2 from formal 
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maintanability démonstratio 


(Vv 


St reports furnished to us by 


tJ 


mremRome Ait Devclopment Center, U.S. Alz Forcs. Mechanical 


Re 


Q 
= 


(Dp 


egtipments included field repaizc tim? pumps used in Fr 
M@etear eClectric power stations (Sets 1-7), and additional 
Memeets Of field repair d2ta from British Aircr2ft, provided 
by the United Kingdom Ministry of Defence, reoresenting 
elapsed times (Sets 21-33), and man hours (Sets 34-46). No 
detailed reports are available for tha mechanical items. 

3. Analysis Approach 


Three approaches were used for testing data. 
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In sets 21-45, the two numbers represent two 
Meme da-a. The first part is ELAPSED TIME and second part 
=s MAN HOURS. 
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(a) Rough check using a technique called Boxplot. 
This gives a rough idea that the data comes 
from symmetric (normal) or skewed distributions 
lognormal, exponential or gamma). 

fou Plotting the data using appropriate probability 
paper. 

[ele Scatistical tests of significance for assuned 
dais >i busvons. 

The following procedure was used for analyzing the 

data, 

i. Sketch boxplot and histogram of data and 
determine if data appears to come from a symmetric or skewed 
Gece tbution. Obtain an idea about center and dispersion of 
data. 

we Oc cle data Of appropriate probabiizty 
Meet for the exponertial and lognormal distributions. In 
some cases, plots were nade for the jyamma distribution. 
322. A Chi-Square Goodness-Of-Fit test was per - 


formec for lognormal and +xponential assumptions, [R¢f. 2] 


56 | 
(D 


eee wi ese, A= tO Shapiro and Melk, {Ref. 3} 
was used to test the lognormai assumption for sample size¢s 


Mess chan 51, (due to availability of the tables), 


WZ 





v. A Kolmogorov-Smirnov (K-S) test was used for 


lognormai and gamma asSumptions [Ref. 4}. 
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II.STATISTICAL CONSIDERATIONS 


pee ROUGH CHECK METHODS 


Purpose: Many batches of data pile up in the middle. 
To analyze the behavior of a batch, we need @ picture of 
where the middle (median) lies and how the tcails relate +o 
Meme ine Middie is generally better defined than ihe cails 
SO we want to see more at the tails. Some values called 
outiiers are so low or so high that they seem to stand apar* 
from the rest of *he batch. 

Ome secinigue £95 displaying such batches of data 
Which is often more convenient than 2 histogram is «he box- 
mee. But ~he histogram 1s still good for showing the shap2 
Meee Gistribution. Boxplot will not show muictipl= modes. 
There may be several "modes", and we want *o ses them. 

A boxplot is obtained by first calculating she lower 
and upper quartiles and the median of the batch (sample) of 
Mameers and then plotting th= numbers cn a horizon«+3:1 lins. 

MWemuewe= GUartil= 1S the value hat divides che 
Meeen L1h<c<O “WO parts, with 1/4 of the numbers below this 
value and 3/4 above it. Th2 upper gqiartile is the value 
Wath 3/4 of the observations below and 1/4 observations 


above it. 
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MiesvexcesteD 2S ©) GradW a NArrOWw rectangular bcex 
with ends corresponding t~> the lower and upper quartiles, 
ama) tO display the median point by 2 plus sign. The length 
Memenis pox is called the interquartile range or H-Spread 
(Figure 1). 

Teo the above are added lines marking 1/2 H-Spread on 
each side of the box. Data values cutside «he lines (out- 


liers) are marked with stars and rectangular circles. 


“Kx tT ly [ ----------- * Ot 0090 


nierquartile rang=2 


a 
«> 2nterquartils range 


interquartile range 
“Figure 1 - Boxplot- 
peast the Gata .ares ordered such «hat, 
POE eels) ss < claiieistceics « S X(N) wheres, 
X(1) : smallest (minimum) value 
X(n) : largest (maximum) value 
median=M , centers of batch 
*-f n is an odd number, M=X((n+1) /2) 


Be on #@S an even numbe> , M=(X ((N/2) +1) +X ((n/2) -1)) 72 


IS 





quartiles=H , center of half of batch 
eeta= (1a (M) t+) /2 where {Al means integer part of A 
interquartile distance = upper quartile -lower quartile 
Here are good boxplot examples for the exponential 


and lognormal model (data generated rcandomly). 


BOXPLOT EXPONENTIAL 


-I + I-------- 0 


=Fagure 2 = Boxplot Exponertial- 
BOXPLOT NORMAL (log is taken then boxpiot) 


ears 3 = Bexrlot Normal- 


2. Histogram (Ref. 35 and 6] 


jv 


Eexpleo: is good fot comparison but a histogram is 
fee good for showing shape. Sometines histogram is more 
Meworul than the boxplot, especially for detection cf more 
amps in the data. Histogram prints data in classes. Exan- 
ples are given below in Figure 4. 


HISTOGRAM EXPONENTIAL 
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HISTOGRAM NORMA 


L ( 
MIDDLE OF NU 
OB 


og is =aken then histogran) 
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BE 
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fee O 1 ae 

ni. 5 3 viel 

2.0 6 HK KK 

2.5 6 tok Kk 

3.0 S tok OK x 

See 2 x* 

4.9 1 * 

aD Z "So 


~Figure 4 ~ Histograns- 


Peo PALTISTICAL DISTRIBUTIONS 
1. The Exponential Distribution 
The probebility density function of the exponential 
Meer spution is [Ref. 2], 


-A 
f£(x,A)=A@e . ; = oO ae 0 


=SeCiMmiLarcive distribution Eunction is, 


~\z_ ~\x 
F(x dy= ‘fhe dt=1-e 


which can be easily evaluated. 


The exponential distribution is used fraquently as a 


{ 
iD 


same to failure model for a system when a constant failu 


ce) 


Beee 2S assumed {Ref. 4]. But for rapair times, it can b 
Shown that it is not an appropriate nodel. Because repair 
Meme inciudes diagnostic, correction, and verification tasks 
arepair time density must have a value of zero at time +=0. 
me Shouid then increase +9 its maximum value rapidly and 


then gradually decrease towards zero as time increases. 
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The maximum likelihood 2stimator of *he parameter of 


mee exponential distribution from given daa is, 


This distribution has many different shapes for non- 
negative variates. Itis skewed to the right, the degree of 
skawness increasing with increasing values of standard déevi- 
ation. The mean and standard deviation are scale and shape 
parameters respectively. 

A random variable is said to have a lognormal dis - 
mmmepit2on if the logarithm of the variable is normally dis- 


1 fi (la y 
F(x)=———— [> 2 CG ax 


<ributed, 


XG 27T 
O 
wlth mean and variance > 
ee. (REX) and OG = Variance(inXx) 


ies da<a i Gs ae Ee Nee #*@e#@e@e#e#e#«@ g en) and aa 
mie neatly =i 2), «..6, L(R)=inxX(n) thar, 
fememam iztk=lihood estimators of the lognormal distribution 


are, 


1S 


Ac 1 n eae 
G=-— 22 (ie jy) 


n 7= 1 
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th 
The p percentiles is calculated by f Refs. 7 and 8] 


Median 2Z(.5)=0 r= ei 

th ; #1.2820 
SO) pers 92(.9)=1.28 2 A : = Rie 

th uti.6459 
Poeper. 24(.95)=1.645 a ae e 


3. Gamma Distribution 
random variable x has a ganuma distribution, and it 


is referred to aS a gamma random variabie, if and only if 


its probability density is given by, [Ref. 9] 


| 

| 1 A= 1 =x 75 

| wren - - -X 2 EOE x 20 
E(x)= a 

| Se Ril) 

; 

{ 0 elsewhere 

where ae> (Oana D2 0 


a is shapé parameter 
b is scale parameter 


Meets fkeongenna funce 107 


iD 


The mean and variance of the gamma distribution ar 


even by, 
mean=a.bd and Variance=a.b 


Special cases of the gamma distribution play an in- 
Memean~ roles in statistics. For instance, for a=1 and b=Q, 


Memget ~he exponential distribution with \ =1/0. 
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For the gamma distribution the scale and shape para- 
meters can be estimated from data in many ways; here are two 


pee Method Of “Moment {M.O2M%) 2stimators 








Z 
(sample mean) Sample variance 
Sample variance Sample n2an 
where mean=a.b and variance=a.b 


be. Maximum Likelihood Methoi (M.L.E.) 2Sstimactors 
This method is more complex as compared to the 
Moms: 22 Lequires numerical calculation. There is an APL 
mumee ion named ESTGAM that gives directly M.L.E. (App.C). 
It is used because it provides a better estimator of ~he pa- 
rameters a and b in large samples, provided «he data «ruely 


come frcm a gamma model. 


eae PROBABILITY PLOTTING 
Ploc~ting data points on probability paper is quite easy 
pmeedees ZO. Legquire many calculations. ft is 2 visuél exa- 
Memeo 2On Cacher then statistical calculation (Figure 5). 
Meares Can we jJearn from a probability plot? 


[ie vestal ind#cact2on Sf asSumed distribuzion. 


2. Examine deparzures from assunecd distribution. 
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gute oa 


Sa. Tails longer than normal 
ee Llaliis SES 5 ee ceca 
-Figur2 5- 


ee LOgnuormal Distribution 


Gonsidqer fhe orisr statistics 


of a sampl 


PU ZS sso oie) ow ow «5 0 e's 0 


Mmemed estimate of the (1/N) th percentile of X 


TE X is N(O,1) then we expect 


-1 
X(i)VHL (i/N) for 


-1 
Then a plot of X(i}) vs @ (i/N) should give 


Seeal.cqh= line. The only problem is 


Bf (eeu) 
PROCEDURES 


a. Use normal probability voaper 


Plot (1/N+¢s) vs X{2) 


Deeeeliccerormolos= in APL “workspace [ Ref. 


Use MINITAB 
me vs@ORe Ss Oo = Cl, put 
meee C2 V Suc) 


NSCORES ar2 (i/N+ts)} 


Zz 


. 
SE Ag, 
- 


—_— « 


c 





< 


| a (er ate | 


-1 
=n niey, so we uss 9 (1/7 (N+s)), 


2 


ea oie 


ZA 4 


eae 





Gomseder oq Of data Y(1) Ss Y(2) §$ 
Mess sss Y(N) where Y(i)=LnX(1) iad then apply above 
Beoceaure to get lognormal ovlot. 
Pee oxpcnential Diszri bution 

Suppose X is exponential A with sample size n, so 
Pr(X < x)=1-exp (-AX) where X > 0. Then (i/nts) is 
an estimate of Pr(x S$ X(i)) =1-2xp (-Ax (i)) 

Semit the data ccally come from an exporential dis- 
mempution (:/N+s) VS: 1-2xp (-AX(i)), should yield a 
eeeeigh+ line. 

PROCEDURES 

as in(i/Nts) vs. X (3) 51 Fectangular paper, 
De oXRONPLOT 2a APL workspace, (— Ref. 27 | 
Ca tO UNTIAE, 
d. Use Gamma prebability paper for ALPHA=1 
BETA=Q where Q is exponential paramet2r so that, 
meen 7 Op ee KO (-k/0) E£Oo> X 2 0 


3- Gamma Distrib 
There is an APL program function named SAM chaz uses 
estimated gamma parameters and gives t25t startistic D. The 


— 


poses. (see Appendix C) 
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ent ro LLCAL FTESTS FOR DISTRIBUTIUNAL ASSUMPTIONS 

We need tO know whether we are actually dealing with a 
Sample from an exponential population or whether the data 
values of random variables have a lognormal or some other 
Meeeeribution. <A statistical test sf significance of a dis- 
tributional assumption provides an objective technique for 
tentatively assessing whether or not an assumed model pro- 
vides an adequate description of the observed data. 

There are three basic steps involved in statistical test 
mexhods. 

fee A tes= Statistic is calculated from the data, 

See rhe probability of obtaining the calculated test 
Statis-ic is determined, 

3. Assessment is made of the adequacy of the assumed 
fest r ibution. 

peeeteneie Dr Obap lL Afy of Oobtat ning ~he calculated 

Meet) Stacistic is "LOW", one can say that th assumed dis - 
tributicn dces not previde an adequate representaticn. 


d 


ib 


ice New pmOus0. 1 1T.y Of Obtaining the calculat 
Meese Statistic "HIGH" , then the data provides no evidence 


that the assumed distribution is not adequate, judged by <zhe 





ipewdecremron a. LOWY or "HIGH" deperds on «he 
user's preferences, and the consequences of rejecting «he 
Seeeerreput. On. For examples, the rejection statistic is tyvi- 


Sey be 0.01, 0.05, or 0.1. For th purpese of this study 


j-- 


> te 
Paes 


ct 
‘@) 


0.05 was selected as the rejec a).p 

memctowid p2 point 2=d Cuz that a Statistical +=st, 
although it allows one to reject an assumption as inad¢quate 
do2s not prove that the assumed distribution is correct. 

ie ihe fee Wess 

This test is ons of the oldest and mest commonly 
used for evaluating distributional assumptions. Basically, 
the given data are grouped ints cells and compared to «he 


expected number of cbserv2tion in the cells b2sed on che 


memamed Gistribution. f[fhen if the test statistic, calcula- 


(D 
a 


Seageitrom this comparison 2xceeds a critical chi-squars valu 
the assumed distribution is rajected. 


iro somey peeplen 25 divedtrg che data inzo cells. 


fmeme-ne Number of observations is small 7 
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#7) 
i or) 
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03) 
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bt 
07) 
or} 
iQ 
Q 
iD 
i 
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ted to use a number of cells as larg= 2s possible, subjec* 


Mmemeene TeSstriction thet :t should not exceed N/5 (N is san- 





The computations involved in the Chi-Squared Gocd - 
ness-Of-Fit test were mad2 by a compiter program, in which 
are given aS inputs the assumed cumulative distribution 
function, the number of observations, «he number of equipro- 
babie cells, and the number of parameters estimated from the 
Soampee. The outputs are the Chi-Square statistic and its 
probability of exceeding 32 chi-squar=> value for a given nun- 


ber of degrees of freedon. 


) 


The procedure used is described as follows: Ref. 2] 

ae The cell boundaries are dérnermined from the 
assumed cumulative distribution as the values such that «he 
estimated probability of the observation values falling 
Within a given class is 1/K. For each class, 


Pe(k=x(2)) = i1/k, where 


iD 


x=the random observation +0 be2 assigned to th 


Mmajeone ach cell Soundiry to be solved irca 
the above formula 
M=nnemm@itocs OL Cells, =f this study k = s/5 


The lower bound of «he fi 


tH 
FF 


Teco sande cne Upper 
meunag Of the last cell are the smallest and «he largest 


values ci the observations. 


ae 





b The expected number of observations for each 
cell E(i) is equal to n/k for each cell. 

c. The number of observed values in ¢ach cell M(i) 
is counted based on the results of the above equation. 


fees the test statistic is 


We) 
iw 


X=- 


+ | 
NO 
| ———j 


a 


e. The computed value x is used to comput2 the 
Mewes, Of Significance, or the probability of a Thi-Square 


value with degrees of freedom equal to K-M-1 where M is «he 


number of estimated parameters from the sample to exc2zed 
calculated chi-square 
2 Z . ; 
A=Pr ( : (Vv) 2 X ) with V degrees of fresdonm 
—X 


If « is less than or equal t> .05 the aSsumed dis- 
tribution can be rejected as inadequate. 
2. The Kolmogorov-Smiznov Test (K-S TPesz) 
K-S test examines the cumulative distripution func- 


Memeeecasher than probability density function . It simply 


- 


oF) 


Sem@pares the sample cumulative distribution function S({.) 


Memmene nypothesized distribution funstion F(.). AS a near 


26 





sure of comparison the test uses, D = Max{F(.)-S(.){3 this 
is the main difference between “he Chi-Squared 5.0.7. test 
muaethe K-S. test. [Ref. 4 ] 

On the other hand here we are dealing with Type I 
Seeger, ~hat is the hypothssis is trues and rejact it. The 
Beepabeia-~y of commiting 2 Type I error is called *+he signi- 
ficance level for the test alpha. Here alpha is *axen to be 
0.05. This means no more than a 5 % chance of rejecting <=hs 
true hyphothesis. 


X has CDF F(x)= Pr(X $< x), 


lA 


Take a sample X(1) S$ X(2) br aetcuewerereler =. XUN 


We want to estimates Pies) = Paeqx < x)= S @& 
n 
s(x), =..9 ESE phew, C1) 
" 1/N "6X > oleae: < ) 
: 2/N i Moz 5 x (2) 7AND.X << 213 
" @ we ® ® ® ® ® ® ® e @ e@ = 
" K/N i Mee hie AND. A  < X.(K +1) 
r : Rs ke ee 
nt 1 ene eth) 
Find D= MAX] S(x)-F (x)| 


Weemg Table 9 of reference 4 For given significances 


Ti D) exceeds the tabulated value that comes from =he 


mame, reject the hypothetical distribution. 
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Pon ono c eo jIGeene caenbulai-sd valu= <hat comes 


mummeca ble, do not tejec= the hypothstical distribution. 


nS 
07) 


There are two Kinds of tables for «he K-S tes-. One 


of these is used if the parameters are estimated from data, 
ed 2f th 


2Te 


Meme ns Ochar one is u 


tn 


a 
— 


4 
42) 
=| 
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—_ 
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oP 


é gjJiven. Here 


(D 
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parameters are estimated from data. 


3. he W test for Lognormal Assumption 
The WwW =£est 1S shown in Ret. 3 £0 be an sffsective 


procedure for ¢valuating *he assumption of normality against 
MennorWel alternatives, evan =f only a relatively small run- 


ber of observations aré availabie. dahn and Shapiro suggesz 


my | 


aaa 2th 


tb 


W test may also b2 used to ztvaluats “h2 assumptions 
meee 2S Oqdncrmal distribution. This follows besaus2e of the 


Beoperty chat if the logarithms of the observations follow a 


wD 


MermMal distribution, then the original values of tha obser - 
Vations are lognormally distributed. 


The following procedur2 was applied using <zables 


a. The cbservatiors are ordered, 
Majer << (2) cess X(n=1) < X (2) 


Poeeine EGLLOW.n 


ve) 
| 
ey) 
ie 
ry) 
= 
(Dp 
Cr 
(WD 
14 
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az2 computed, 





ee B 2 
hehe (ae 
2 2 Z fi=7 (-) 
Ss 2, [12x] - 





eto sws Ven, k=n/ 2 
pees Saag, “K= (n- 1) 72 
k 
B= (meek te yl nx (n-2+ 1) )— (inx (2)) 3] 


oe 
Weese on] Values of a(k=i+1) f£5r i=1,2,..k 


are given in table IX (Ref. 2] for n=3 to n=50. 
Mw mbHe CoS Sta-aStic Ww -is 
2 Z 
N= 9527 5 
Gee ihe a2pproximac] probability of obtaining «he 


[mmeata-ed Value of W can be obtains? from the formula 


|w-é | 
z=Aspin |--=| 
,i-w| 
using the values of A iceceg= ven ine Fable XE of Ref. 2 
for the appropriate? sampl2 size and then using standardized 
Gemmal distribution to determine the probability of obtain- 


@ygea Value less than or 2gual to Z, which is the sigqn:ifi- 


cance level of the test. 


G@e= Pr{Z Ss Z} 
gq 
femg@uis less than of equal => .05, the select=d ie- 


femme Sign2-ficarce in this study, the lognormal is rejected 


as an inadequate assumption. 
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1. Summary of results 

Tabies 2,3 and 4 summarize the results of the sta - 
tistical test analysis. In all cases the 5 % level was che- 
Pemeas the reject criterion. The 46 sets cf data used in 
this study can be divided into three parts as indicated in 
Section on IB2. 

The following is an example of the computer procram 
Summary table which ircludes the results for the Chi-Squares 
test and caiculated parameters from +the sampie data. 


SET NO 18 - AN/GXS-2(V) 


nab apwe Eee? Ee ee oP Oe eee eee Pee Pe eee eee eee eee eee eee eee ee Pee ee eee eee eee Se ee ee ee 


SAMPLE SIZE N= 26 NWOee OF *GSELS ~ kK = 5 {a) 
SAMPLE MEAN = 19.70 STANDARD DEV = 22.62 (b) 

EXPONENTIAL  LOSNORMAL = DIFFERENCS 
DARAM1 0.05 2.56 (c) 
PARAM2 ORR Sa, 

i ea hen 0 Ss 2,96 +% (cd) 
PU-=tH PERCNT 56:5 Zao 3 = elo gee bm 
Both PERCNT oo eae ae 4Q.22 l2.7 7 
Pes tH PERCNT 59.0090 55.46 Gao Ve 
CHI-SQR STA? 5.92 0.54 (<) 
DEG OF FREED 3 2 (f) 
eau. F LEVEL O29115 E+ C00 on Ou O (gq) 


Dap @P <P aP «ep <P <P ap mana 4 Pa a Se ee 
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TOTAL DATA ANALYSIS RESULTS 


GAMMA 
eS © 


ix 


Gn. 
Sa. 


NORMAL 


LOG 

eo 
BO. Fr. 
Pil. 


ick ose) 

Ana Syae 
| aaa OF 0a 
Sd. 
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Sample 
oy Z| 


S= _ 
No 
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Poles 
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! 


05 AS 1 A 0 GG TE te tT 0 0 et 0 0 tt 0 0 0G ot et GG Ot 0G tt att OG 


+H+H+ H HH + ++ 


OEE (OG COIN MOG OM NOVO DMM eRe MOOG KO NM eT MOM OG Ne OOF MG ei Cri tet Tet Cosa A 


Se NR ee eee He ois i ale: (4 ytet fs fii a fue 
1 I 


Mn a ase ood Oep een eet eect. Pest ect ort i set ork i padeat—ieetess eel *<Peek ok ko aet otek het aol ot eo i oor as 
OOCOMR MMM Oma MO BM use COMI9OMO MAb MoMO BOM MmoOocrsoorsmrsm 


QOOUMM MAW OmAUAMURMAUOUNMOUMAUUMMUMVUUAMUMmMUUUmAOM 


SHEE HE HHH HEHEHE HEHE HEHEHE FEE FHEHEHEHEHEHE HEHEHE HH FHtte + FH tt 


Mmm Ae OM OMe mato Meet OO mm OM Me tM mG et OG ANG 


MAMMMMM NWA. AMMA, OM 0) i M M Ay fey OO fa hy OAM OM OM fy My OO Oi fy Fry MO Fry fa 


MMOAMNMONK HODHDOMODOMOOOHF FH DODONSKMEMNK AMF KK—OOMNWN OMrnr 
FIMNANAYIA AT NNT MAMMN ANIMA MDANAD ACN ZF MOASTMH MMO AGHAMLIAM ON =F 
Spm 


w 


ANN FINO DHNOK ANN ANUOP DHOOM AMs+ 
NN QION CIN 


o 
~~ @ 


less than .05 
bulated valu 
beg OL 2Ox DLO. 
Good 


Pr (WwW) 
G 


eon 16 cl 


g 
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LOGNORMAL AND EXPONENTIAL RESULTS 
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Tadpulated 
Value 


A/R 


st) 
value of 
ws. ES cs. 


= 
= 
— 


D 
M.L 


A/R 


Vaiue of 
wW.O.M.ESt. 


D 
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ae The number of 2quiprobable cells, K, was chcsen 
as N/5 where N the sample size. 

be The sample mean and staniard deviation are cal- 
culated based on the maximum likelihood estimates. 

eee Panam ly for the exponential distribution is ths 
reciprocal of the sample naan, for th2 lognormal distribu- 
mon PARAM! and PARAM2 ars dD and GS Pies pad=allece ns  O= he 
megmorMal distribution. 

d. The MTTR and the 50th, 9)th and 95th percentiles 


are caiculacted from the sample and are based on the rela- 


4. 
O 


—_ 
Nae 


nships between the calculated parameters and their dis- 


Ce 
ty 
| 


ibution functions. The percentage difference is between 


ct 
om 
D 


= exponential end lognermal ATTR and percensilés. 
pee he Chi-Square statistic. 
a Che nimber Of degrees of Freedom (K-M-1) is K-2 


Meemcone exponential distribution and K-3 for «he lognormal. 


(p 


This is because one parameter is estimated from “he sample 
in the =xponential case and two in the lognormal case. 


Spee te wee Velh Of Sslgnlitcanc= 2s “he probabiiity of 


@ Chi-Squar2] variate with the specified degrees of [reedor 


Smeecding the calculated Thi-square sta+istic. 
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Po =n ex DOncicitt model, Only the Chi-Squered test 


Oeeece-Ue Gane exPDOren- 


f- 


was used. In 6 out of 46 cases we wou 
tial model. If w2 compare the exponential vs. lognormal no- 
del using the Chi-Sgquare test, test orefers the lognormal 
model tc th2 exponential nodel. In only one case would the 
exponential model be accepted and th2 lognormal rejec-ed. 


exponential 
A R 


A 5 19 | 24 
legnormal Sa | =| 
FE ; 1 | 21 | 22 
6 40 


Therefore the Isgnormal model is a becter descriptor 
With respect to the exponential model. 

In most cases the lognormal 2ssumption works reason- 
memye well. in those cases in which the results of the chi- 
Smee d teS=, W ~es= and K-S *est do not agree, a plo«= on 
lognermal probability paper was used to determine the appro- 


priateness of this assumption (* in Table 2). Comparison 


(Db 
(Dp 


between *est methods were nade in thr ways for the lognor- 
mal model and cabled below. 


@e Chi-Sq. vs. K-~S test (46 cases) 


K=5 
A R 
: A pots | 10 4 29 
Chi-Sq. |--=- |---| 
R | Sid || 21 
18 28 
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be Chi-Sq. vs. -W cest 
a 
A R 
ae Vee | 21 
_" tf &© | 94| 13 
23 11 
Cs Wvs. K-S test (34 cases) 
Kae 
A R 
A a3 10 | 23 
W eels ot 
R i 3 | 8 | 11 
16 18 


For <th2 lognormal methods the 


give different results. 


test agres 
the W test 
Therefore, 


three test 


(34 cis¢s) 


£ mehods sometines 


-- 
Se 


The K-S =esst and the chi-square 
in 74 % of the cases. The chi-square agrees with 
in 82 % of the sases. @ and K-S test agres 62 %. 
there is no significant difference anong “he 
methods. 
A previous study has been dene using 24 data sets of 


electronic systems and equipments, using only the W and chi- 


squared tests for the lognormal model [Ref. 1]. The cabis 
below gives a comparison between the 2xponential and lognor- 
Secor che previous study. For the axponential model do 
not reject (A) or reject (R) decision was made oased on the 
chi-squared test. For the lognormal model accept ay either 
the chi-square or W test was used as 2 criterion. 
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Semen. +82 


A eae 1h | a 
lognornal | “5a tae 
R lie Zl 3 
8 16 


The tables below show same kind of results fo 


ty 
ct 
+ 
iD 


current study. Only diffsrence between the two studi 


(ID 
~” 
j 4 
WV 


there are different kind of systems and equipments 


ib 

6 
'Q 

& 


mechanic or electronic) in the current study. 


Genco eset 


A l-> | 22m 28 
Logner mal Pa anes 
R 1 fle 18 
6 40 
Wo creed se Seno Cases} 


vo a 1 8 
Lognormal {----|----| 
R | 90 ¢ 351 
1 eZ 


For electronic 2quioments only (13 casés) 


These results indicate that the lodqnormal model is 
better than the exponential model for repair time daa. 
For the gamma model, she MLE estimator is more po - 


Metrul than the MOM estimator, with the K-S test used for 


testing the data. The dacision was mace *o acce 


i 
ct 
vt 
aay 
(p 
t4 
iD 

( 


Mees Of che MLE when the two estimators gave dirferen* 


ay 


s 





results, since in no case did the MOM suggest acceptance 
when the MLE said to reject. In the sight cases where thers 
was a difference, therefor2, the MLE was used as the crite - 


rion. This gives the following results for the gamma nodel. 


td 


MOM 
A R 
Pear | 8 |{ ZY 
MLE j--=-|----| 
Pee Ost A Sarl 19 
19 27 


Therefore using the MLE, we would accept the gamma 


rt fa 


mem@ere an 2/7 out of 46 casas. Ir 16 of these 27 cases, *he 
lognormal model is aise accepted and in 5 additional cases 
all three distributicn assumptions would be accepted. Again 
only one case (case 6) would we accept the gamma and expo- 
nentiai model but reject the lognorm:l. 


We can also compare lognormal and gamma modéi using 


mae K-S test. 


gamma 
en Ee 
lognormal ----|---- 
fies to 28 
26 20 


The above table shows that the K-S test prefers +he 
Gamma nedel. 27 out of 46 cases accepted the gamma model, 


Dut only 18 out of 46 cases accepted the lognormal model. 
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Tables 6-9, (Appendix B) Sumnarize the results for 


all cases and for each of the three different types of data. 


t 


Tabie 10 (Appendix B) gives total results in terms cf sampla 
eeewes Of 20-29, 30-39, 40-50 and over 50. For exponential 
assumption, increasing sample siz2 increases rejection per - 
centages. For gamma assumption, increasing sample size also 
increases rejection. In lognormal, increasing *h2 sample 
size decreases rejection up to sampl2 sizes of forthy, “hen 
rejection increases rapidly. 

Because of the large number sf data sets, only some 
interesting cases of the data analysis and probabiiity plots 
for those sets are discussed in the next section. 


2. Discussio 


ts 
{O 
(th 


specific ¢ 


is 
(tA 
iv 
yn 


The following discussions illustrate some cases 


th 


Which are of special interest. The data analysis results 
and plots are given in Appendix A. 


Pe oSe NOs ero — Ba GXS~ 2 (V) 


This case is én illustration where all three 


im 


models azr¢ accepted (not rejected) for all tests. Data wers 
memlected from 26 maintenance tasks. In this case, 4s in 


most other of «he formal maintainabiliity demonstration {(ca- 


ses 8-20), the tests were performed in accordance with the 
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U.S. military standard for maintainability demonstration 
(MIL-STD-471) [Ref. 23] in which sample selection was ge - 
nerated from knowledge or astimates of component faiiure 


rates and tests were conducted uSing service *echnicians of 


ct 
is 
iD 


the appropriate skill levels and given prior training on 
equipments. 

This case fits the lognormal, expenential and 
also gamma for both estimators. All three tests give high 
pagmaticance for lognormalitv. 

pape ses NO. © = Reactor Coolant Pump 
This case is the only one the exponential model 


accepted and the lognormal model rejected. However th 


wW 


plots are not good. Further investigation is necessary in 
this case but since only the data was available, it was not 
Peestbie for ws +o de so. 


ce. Sets No. 28 and 41 - Main Rotor Head and Biliad 


(b 


S 
These two cas2s indicate a very good example of 


Memnermality. All three + 


wD 


eS Siena “dorwnot rejec:" fcr 
the legnormal assumption. The plots are also very good. 


sent «zhe data sets. 


(D 


The gamma family does not well repr 
dj Sets No. 7 and 2 - French Fessenheim Pumps 


( Down / Repair ) times 
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These two cas2s aze discissed as unsuccessful 
for gamma and exponential but good for lognormal since art 


least 2 out of 3 tests give acceptance for the lognermal 


iD 


assumption. For repair times the W test rejects th2 assump- 
tion while the Chi-squared and K-S tests do not, but down 
imes are rejected by only the K-S test. There are some 
Outliers around 5%, but the rest of the data show a gocd 
Meagnormal £22 in the probability plot. 

The existence of outliers could be a reason for 
Mmemection by one of the 3 tests. For downtime data, after 


Trimming the two high outliers the lognormal assumption was 


Soe Ourreibaic Cie, a@iee> trin- 


()) 


also accepted by the K-S ¢ 
meg One Outlier the data shows high significance (0.15) for 
the W test. Thus, in som= cases it would be valuable «+o 

€Xamine in more detail th reason for «he high or low outli- 
ers. In this case, only the data was made avaiiablée so that 


emas Was rot possible. 


Sees eNOweed) —- Haater Drain. Pump 
Set No. 5 - Main Feedwatesr Pump 
Set No. 4 illustrates the case wher2 all thres 


tests for the lognormal model are rajected. Here outliers 


az2 important factors. The effect of outliers can be seen 





femely from =he boxplot, nistogram and probability plot. No 
5 is almost the same as no 4, the only difference between 
the +wo being that the K-S test barely accepts the lognor - 
Meaty in set no 5. For these cases, “rimming the outliers 
did not change anything in so far as significance levels or 
acceptance or rejection are concernsi. 

f. Set No. 9 - AN/ARC-162 (Intermediate) 

Set No. 10 - AN/ARC-163 (Intermediate) 

These two cas2s are discussed together, because 
they give opposite results. Detailed reports ¢xist in both 
cases. Test technicians were used for running the tests and 
removina the faults. [Refs. 12 and 13} 

The exponential assumption is rejected in both 
cases. In set No 9 the lognormal assumption is rejected only 
Myeene K-S test while the other tests give high signif:i- 
Memce. The probability plot and boxplot also signal «hat 
Sao not reject" is appropriate fom this case. The gamma 
Meme! fits this set for both estimators. 

MawoceeN@. 10 only the W test accespis, bux the 
Chi-Squered zest barely rejects (0.045) and the K-S t¢st 
rejects the lognormal model. However the boxplot and proba- 


meet y plot are fairly good for this sez. An 2xaminaticn of 
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the report indicates that some bias 2xists by the use of 
"average" times for some r2pair time segments for much of 
mae date instead of normal data. Therefore, a fusther exa - 
Maemeation of this data would be in order. 

ge. Set Now 15 - AN/UYK-14(C) 

Set No. 15a - Revised AN/UYK-14 (C) 

Set No. 15 illustrates what happens when the da- 
ta points are not selected at random and thus introduce bi - 
Meseeern this case, 33 tests were selected according to the 
Sample selection procedure in MIL-STD-471A [Ref. 23]. How- 
ever in order to have 50 data points according to the test 
plan 17 additional data points using "average" repair «imées 
for several components from a previous test were used. This 
therefore introduces a bias in the data, resulting in the 


mee lal rc 


uw 


Necee onmofesth= 1oqnormal nodel (Table 3). after 
taking out <he biased data points th? analysis shows that 
mies LeCqnormal model =s a good descriptcr. The histogran 
does not look very good, but the probability does for the 
lognormal model. 

This exampl= then shows that if data does not 


Mee an assumed distribution and if the reason is bias, then 


a | 


@tter finding a good reason for throwing out some bias po - 


& 


Pres, the data frequently will agree with a lognormal model. 
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he Sets No 21 and 34 - Fuselage Structure 
These two cas25 are example of large sample si - 
mes) (157 and 171). For both data sets all distributional 
assumptions were rejected by all tast methods. But for the 
MegmoOrmal case probability plots, boxpiots and histograms 
are fairly good. These cases should be further investiga - 


ted. 


Pee kkOKS IN CALCULATED PARAMETERS 
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The percent error in the maan-*time-to-repair (MTTR) 


meecealculeted as follows, 


etec ees.) 

O¢ ax 

E = J ae ae 109 wheres, 
4 
Lo ~ “ 

+0.56¢ 

M = the lognormal mean = 3 

loa 

M = the exponential mean = xX 

exp 


The results of ths percentage error in the mean 
which are summarized in Table 5. All cases neve an error 


Mss than 18 %. Leenoseso Somes 2£he =nheren 
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formula normally used is, 
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mmamsaince £0©r a reasonable availability MTTR << MTBF (e.g. 
MTTR is measured in hours, MTBF is 100's of hours), then it 
is shown that the exponential assumption for MTTR introduces 


negligible error in availability. 


eee CONCLUSIONS 
From the data analysis in this study as well as the pre- 
meetemotudy { Ref. 1], it is concluded that the lognornal 


model is a good descriptor for repair times. 28 out of 46 


pb 


Sets show that, within 2a 0.05 level of significance for 
least one of the 3 test methods , this model can not he re- 
jected. Similarly, the data analysis shows that the assump- 
meee or tne exponential distribution should be rejected in 


40 sets. The gamma assumption is not fr 


w 


jected in 27 cases. 

The percent differenc>= in the MTIR, when assuming ¢xpo- 
Menetal distribution when the true distribution is licgnormal 
has been found +o be negligible for calculated system avai- 
lability. 

Mmer the results of the statistical tests indicate oppo- 
Mmeemconclusions, probability picts, boxplots and histograns 
are helpful in determining an acceptable model and =o explo- 


re anomalies in the data. Prebability plots are also useful 


we) 
2 
wm 


stimating percentiles, as well as density parameters. 
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Seow iaa= REVISED AN/OYK-714 (C) 


SeUPLe SiZe N= 3 NOwsoF CELLS ~xk = 5 
SAMPLE MEAN = 14.37 STANDARD DEV = i 66 
EXPONENTIAL © LOGNORMAL = DIFPERANC2 
PARAM go. 07 2 
PARAM2 Ot 
oA ARy 14.37 14.19 Vere os 

ego rH PERCNT She Sits 5) ares hes, 7.386 
39Q0-TH PERCNT Bro 18, 24.83 63450 % 
Pomrd PERCNT 43.06 30.40 41.65 % 
Si-SOR STAT Boni Pog es: 
DEG OF FREED + 3 
SeGNIY LEVEL O22 16 E- O04 3.484 2+00 


BOXPLOT SETI5A 


BOXPLOT LUNSETI5A 


-----I + [--------------------- 
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fo UOGRAM SETISA 


MIDDLE OF NUMBER OF 


INTERVAL OBSERVATIONS 
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Times (nia) 
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Set No 45a 
Lognormal Plot 
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SET NO 18 - AN/GXS-2 (V) 


SAMPLE SIZE N= 26 NO. OF CELLS K = 5 
SAMPLE MEAN = 19.70 STANDARD DEY = 22.62 
pee UBC EORHAL See eae 
PARAM 0.05 2.56 
PARAM2 0.78 
MTTR 19.70 1c We 2.96 4% 

50-TH PERCNT 13.65 12.93 5.61 % 
90-TH PERCNT 45.35 40.22 7 TAS, 
95-TH PERCNT 59.00 55.46 6.39 % 
CHI-SQR STAT 5.92 0.54 
DEG OF FREED 3 2 
SIGNIP LEVEL 0.115 £+00 0.764E+00 


Soxe LOL SET18 


--Il + [------ 0 0 


BOXPLOT LNS£ET18 
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HISTOGRAM SET18 


BeDDLE OF NOMBER OF 
EH PAL Gee Tas 
10. 14 Se a ee KKK 
Ze 5 KK KKK 
210) 2 dere 
40. 0 
506 0 
60. 1 * 
70. 9) 
80. 0 
90. 2 *% 
HISTOGRAM LNSET18 
ae OF NUMBER OF 
PATER VAL OBSERVATIONS 
1.0 1 * 
ie 3 er 
pO 6 HK KKK 
ye 5 6 KKK 
B20 5 fe tok 
Sips, z hag 
4.0 1 * 
4.5 2 x 
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_L ognoemal Plot ° ° 
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Seema 2! = FUSELAGE STRUCTURS 


SAMPLE SIZE N = 157 NGS OF CBELS K = 2% 
SAMPLE MEAN = 1.54 STANDARD DEY = ees 
EX PONENT IAL LOGNORMAL DIFFERENCE 
PARAM! els) =O. 06 
PARAM2 Oo 
Del a 54 1253 1.01 & 
Bu-tH PERCNT tOe7 0.94 Tas Si Zelt 
90-TH PERCNT S20 3230 6.93 % 
pain PERCNT 4.63 Pa os, 2.64 % 
Sae-oOR STAT iS )5 Os 55 e) 
DEG OF FREED ce aa 
BeGhiIF LEVEL 0.0 0.0 
Beee LOT Sit21 
-=f 6 #ecsae a 5 OS Or Go 6) CO 
Bere LOT LNSST21 


~--- -- ---- -- ------] + 3E 
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Hes tOGRAM SET 21 
Pac * REPRESENTS 2 OBSERVATIONS 
MIDDLE OF NUMBER OF 
INTERVAL OBSERVATIONS 
om 33 ek ede eK eo OK KK 
1. 73 ee $C ee aK eK I i oe a a eo EK OK 
De Dal eae HK ek kK KK 
3. 7 KK KK 
ty 4 xx 
Je Z - 
6. 5 KK 
yi 3 KK 
8. 1 * 
Qo. 0 
10. 2 % 
HISTOGRAM LNSET21 
MIDDLE OF NUMBER OF 
INTERVAL OBSERVATIONS 
= 3 oceale 
2.0 2 sie 
-1.5 2 KK HK HK KKK 
-1.0 16 Hk OK KK Ok KK KK 
-0.5 34 He He 2 oo a oo i Kok KK KE 
9.0 36 ee HAC 2 I Ck aK ae KK KK KK RK KK 
0.5 27 HK HK eH KK KK KK 
1.0 10 eK KK 
1.5 7 eK KK KK 
220 g KKK KKK 
209 2 *z 
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Saaevgezo = BATN xOTOR HEAD AND BLADES (ELAPSE TIME) 


SAMPLE Sze N= 33 NOweonr ChELsS: K = 6 
SAMPLE MEAN = ee OT STANDARD DEY = a se 
EXPONENTIAL © LOGNORAL ©» DIFFERENCE 
PARAM1 0.48 9.46 
PARAM2 0 Gee 
ead Tames es Jaen bs) 0-596 

s0= TH PERKCNT 1.44 12,58 SAO ome 
EO-TH PERCNT ieee a 4.00 tes ea 
ea-TH PERCNT Bec | DAY i533. 
Sa =5OR STAT ale. o 4 GS, 
DEG OF FREED 4 3 
Beane: LEVEL O70 2-01 Jeo o4EF+09 


Eee LOT SETZS 


--I + I---------- 0 


Beee LOT LNSET28 
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Ste Omste— ur USELAGE STRUCTURE 


nmap 4p 4p ape 4? 42a awe ee 42 eee ee ee eee eee ee | eee eee ee ee eee eee eee eee eee ee eee eee ee 


SAMPLE SIZE N = 171 NO. OF CELLS K = 34 
SAMPLE MEAN = 2.83 STANDARD DEV = 4.60 
eee | LOSNT RE eee 
PARAM 0.35 0. 36 
PARAM2 15 WE 
MTTR 2.83 2.68 5.26 % 
50-TH PERCNT 1.96 1.43 36.77 % 
90-TH PERCNT 6.51 6.03 Meehee 
95-TH PERCNT 8.46 9.05 6.51 % 
CHI-SQR STAT 126.65 92.85 
DEG OF FREED 32 37 
SIGNIF LEVEL 0.596 E-07 0.596E-07 


BOXPLOT SET34 


18 Os oak ae) 018, OIC 9 0 0 


BOSPLOT LNSET34 
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HISTOGRAM SET34 
BOB SERVATIONS ARE BELOW THE FPIRSI CLASS 


MIDDLE OF NUMBER OF 
INTERVAL OBSERVATIONS 
0.500 38 Seok ea ok eo kK ek eK KR EK RK KK KEK KKK 
1.000 33 $k kk feck tk eK KKK KKK KKK KKK KKK 
1.500 25 KKK RK KE KKK ERK KK KR KK 
2.000 15 Ie RHO KK RK KK 
2.500 9 eK Ke 
3.000 13 eK KKK RK KK 
3.5G0 2 i 
4.000 3 ae 
4.500 3 pees 
5- 000 3 == 
5. U0 ) 
6. 000 5 EEZEE 
be 00 1 * 
7.000 1 * 
HISTOGRAM LNSET34 
MIDDLE OF NUMBER OF 
Tene R OBSERVATIONS 
~ 2. x 
eX 
KK EK * 
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Soot tea MALTN ROTOR HEAD AND BLADES (NAN HOURS) 


SAMPLE SIZE WN = 35 NG= OF .CrLis  K = 
SAMPLE MEAN = Be 50 STANDARD DreV = gus 4 Z 
EXPONENTIAL © LOGNORMAL = DIFFERENC2 
PARAM1 si) 2S, OSG, 
PARAM2 Om t 
col gg Be a8 SD) Ope 0a 
Se-tH PERCNT DRI SS Zot We? a 
po-2H PERCNT 550 0 ee liiere Ou 6 
yoean PERCNT 19.49 SIS 8 Gee % 
Sie Ok STAT | US PAY 
Who Or FRED 5 4 
SIGNI? LEVEL = 9.476801 9. 1262409 
BOXPLOT SET41 
——— : : 


BOXPLOT LNSET41 


-~-------- I ’ I --------------------- 
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HISTOGRAM SET41 


MeeDD LE OF NUD en OF 
INTERVAL OBSERVATIONS 
o. 5 4K He 
D 13 tek kk ok ok 
4 10 etek tok ek 
6. 4 tek 40% 
8. 1 “3 
GE 1 * 
le 0 
14, Q 
16. 0 
18. 1 i 
HISTOGRAM LNSET41 
Peover OF NUMBER OF 
PN TERVAL OBSERVATIONS 
— (es 3 jc 
0.0 7 wk KK 4k 
0.5 7 fee a KK 
1.0 3 eS 
1.5 9 tek feo a kok ek 
D0) 5 sek ek 
Bes 0 
3.0 1 * 
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APPENDIX B 


: 


Lar 


ble 6 


TOTAL RESULTS (46 cases) 


_ 22a eae ae ee ee eee ee ee eee ee ee eee eee ee 


BPepeonen=.al Distriblition : Only Chi Square test used. 


Accept Reject 


m_apP 2a ae ane eee ee eee eee ee eee 


{ 
i 6 40 
| | 


= =P ee See eee eee eee ee Se eee Se 


Lognormal Distribution : Chi Square, Kolmogorov Smirnov 


and W test used. 


_ i= 0 No W *est 
A R 
a test wW test 
A R A R K-S test 
| is a ine A R 
A 12 1 7 1 emer erence ee- 
Sie od . | | Ch. Sd. ot Z | 2 | 
test ee ie oO 1s | 
Ray 1 2 3 74 eee ene eae 
{ | 
( 34 cases ) { 12 cases ) 


Wee. =Ss.LnAaOoLls 


A R 
| | 
M.L.&. A | tT? 8 | 


— { 
estimators R j 9 | 19 
{ 
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PUMP DATA RESULTS (7 s2¢s) 


Exponential Distribution :; Only Chi Square test used. 


fearcrmal Distribution : Chi 


5 
and W test used. 


ag a aa rl NOOO wy 2]eS= 
A R 
w vest W test 
A R A R K-S test 
| ee A R 
; A | 0 1 2: Of, eee nererece- 
ene Sq . | | Cie Sa 6) of 0 0 | 
were ee ee eee eo eee test amen | ---- 
i> © i | R 0 1 | 
Ref 0 1 0) 2 (eee renner e= 
Jo SS ee | 
( 6 cases ) ( 1 case ) 


Gamma Distribution ; Kolmogorov Smirnov *+sst used. 


MeO te = Sema ors 


A R 
i 
Mee 2. A of 1 
| | 
| 
est2mators & | 0 | 6 i 


a> 





_ = ap oe ae a ee ee ee eee ee eee ee ee ee Oe ee eee eee a Se 


mam =e ea ae ee ee ee ee ee ee eee 4 oe ee ea oe 


imeqnormal Distribution :; Chi Square, Kolmogorov Smirnov 


and W test used. 


K-S *2est 


ee No W test 
A R 
Ww test Ce meio 
A R A R Ke Sos 
= sea oo SO A R 
| | 
A 4 0 a 0) a ri 5 ge 
Gi. Sq. { | Clee Sa QO | | 
a teen nn rn ee --- test | -----|----| 
test { | R | {i 14] 
R | 0 0 2 Ce | ee 
{ 11 cases ) { 2 cases ) 


Gamma Distribution : Kolmogorov Sniznov test used. 


MO... 2S*tinatcrs 


A R 
{ | 
i ot lai) he A { fi 3 { 
eee eal 
| | [ 
estimators R { 9 5 : 
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PReeO DATE RESULTS (26 secs) 


—e ep ep SP ee eee ee ee ee ee eee ee 


Beepcicn+tial Distribution +: Only Chi Square test used. 


Accept Regese= 
| 
4 i 22 


iegnormeal Distribution :; Chi Squares, Kolmogorov Smirnov 


and Wecos= wesc, 


K-S test 


oo aa ae No W test 
W test W test 
A R A R hao ease 
~ -=---- ------=------- A R 
{ | 
| | eee rn 
Gat Og. { Ciesgs. sAtme ce | 1. { 
- ------ ------ -------- ‘Ceere |----- ----| 
test { { R | 0 | 6 | 
= || ee 1 an poe 2 
anes een 
eases.) ( 9 cases ) 


Gamma Distribution : Kolmogorov Smirnov <tsst used. 


MeO... 2S5S*1imators 


A R 
| 1 
Me ee Ss A { 11 H 5 
EE er 
: 
estimators R 9 10 i 
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TOTAL & 


For exponential dis 


Sample Size 


Z0= 29 
30-39 
490-50 
GT. 50 


mes. tognormal distr 


Sample 
Size 


20-29 
SOs) 
40-50 
GT. 50 


For gamma distribut 
Sample 
Se.7 = 


PAO ae, 
B0=39 
40-50 
GT. 50 


Palermo 


nou L TSe(W.r. 


Serotec Oon  ° 


S53 


pe Ons = 
Chi-Sq.test 
A ik 


53,3 46.2 
SiGe) 335.45 
55.3 os 7 
256 Jase 
aloha: 8 
"pst 
84.6 
Sie ee 
25s 
NG? 


25 


+. sample sizes) 


test K=S2e2 = 
R A R 
Pee es Oe2) 095.5 
Ver ee SO. S10) 
Ul) sooo. | 
--- (S27 S36 3 
Avin MLE 
R in MOM 
ore | 
163-7 
125 
66.6 





ABPENDIX C 


COAPUTEA =e nCOGRANS FOR DATA ANALYSIS 
1. LGNRML (FORTRAN) [Ref. 1] 

The program makes a Chi-Squared Goodness-Of-Fit test for 
the exponential and lognormal assumptions. The main program 
computes expected values of ordered observations, sample 
mean and standard deviation based on che maximum likelihocd 
estimates of the varameters of the lognormel distributicn, 
Calculates the 50th, 90th, 95th percentiles of the exnonen- 

e 


tial and lognormal distribution assumptions, ver 


rh 

ct 

a 
p 


differences petween che calculated parameters o 

Sexponen=ial and loqnormal distributions. 

2. BOXPLOT (MINITAB) [Ref. 5] 
Semcax; B30xDLot Vector 
Description: Generates a boxplot display for a vector 

Semeea, A crectangulaz box with ends corresponding to lower 

and upper quartiles is presented with the median marked with 

an plus sign, and lines are drawn on each side of the bex 
with crosses marking the lowest and highest data values wit- 

Peegeeicdes Guattile distancSs in each side. Outliers are mark- 

ed with stars and circles. 

3. NORMPLOT (APL) [Ref. 21] 


Syarax: Normplot vector 


Pameme Sts: Wid- controls the horizontal size of display 
Dep- controls the vertical size of cisplay 


Saubprogranss NSCORES, PLOT 


Deses.DtLON: Genera =2=S a normal orobabilaty plet for a 
Meeeeot ©f daza. [f data fits normal then p 
bine. 
PeeecoXPONPLOT (APL) [ Ref. 21} 
Peeeax:  Exponplot vector 
PeGmameeters?: Wid= conrtzrois he horizontal size of display 
Dep- cone rolls the vertical size of 4 


Ssmbprogtam: PLOT 








Description: Generat2?s @ exponential probability ploz 
Momma vec OL Of Cata. If data fits exponential then plot 
femoeralight line. 
oeeeeGaAud (ADL) 

Syntax: GAM vector 

Subprograms: ESTGAM, GAMMA 

Description: Perform Kolmogorov-Smirnov Goodness~-of-Fit 
Beatmania Gives test statistic and "reject" / "do not reject" 

a 


for gamma assumption fcr beth MOM and MLE sstimactors. 


Ov 
r 


KOL (APL) 
Synzax: Kol vector 
Subpregram: MOM, INVV 
Description: Performs a Kolmogorov Smirnov Goodnees-OFf- 
Pat cest for a data vector for lognormal assumption and 
Beemer c]est Statistic. If test statistic is bigger than the 
tebulated value "reject" the assumed distribution, otherwise 
Magmmot reject". 
Te W (APL) 

Syocaxs A W Vector 

Description: Performs a W test for the lognormal as - 
Slee on. A must be a data vector (from Table IX, Ref. 2) 
and X vector that is testing against lLognormality. A and X 
must have the same dimension, cutput will be W test statis- 
tic. Then using the value of A,juani € (Table XI, Ref. 2) 
for the appropriate sample size and than using standardized 
normal distribution determine the probability of obtaining 
value less than or egual Z, which is the significance level 


Of =he test. 
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eSretn? 


9L0 
LijsFke| (J-Kk) 

Rie eC pry 

' TEST ST&TISTIC FOR MOM ESTIMATORS 
SLSxX.CRRELI (=) 


9LS&, 


TIAGO NM ANMMM MIAMI HAMNANAMAMANMMMMMAMMM 
tO Cd ON CO GEA OILS POOP SRO PO Po PoP Pe Re RRP ee ee OO NO CS OP 
ON RUS Oe Oe NG Op Gir OOOO Ol b ir Ghibli ht 


| oe ee ee ee es es ee es ee ee ee ee ee ee ee ee ee ee ae ee ee ee ee 


eked 
? 


GAM SET{ 


L7?? DO HOT REJECT FOR ALFHAXO,9O5 FGF 


QA FROGRAM FERFORMS KOLMOGOROV SMIRKHOF TEST FOR: 
QA GAMMA ASSUMPTIORM FOR MOM GHIO MLE ESTIMATORS 


IS ',eRROCL] 


MLE 


foot stHt LT STic FOR MOM ESTIMATORS £5 9.363 


REUJECT FOR G@LFHA=0,.05 FOR MOM ESTIMATORS 
TEST STATISTIC FOR MLE ESTIMATOR IS 9,154 
REJECT FOR ALFHAXQ,05 FOF MLE ESTIMATORS 


? 


( REJECT FOR ALFHA=H=0,05 FOR MOM ESTIMATORS ! 


L3z' 2O NOT REVECT FOR ALFHAXO,05 FOR WOM ESTIMATORS 
wey; irteri+i 

SL GX 1 CTT >t) 

AD? GAMMA HOLIIT 

JSJIe dd y F 

als 
Bee etc is—* ) 

S5¢SC yo] 

Peo estat LS>TiC FOR MLE ESTIMATOR i5 ‘'s,455[1] 
gas et 1 ak= ) 

Resear oe ALehiaxO, OS FOR MLE ESTIMATORS '* 
Do? 


EoaTrixatlors 





9 KOL COI 
7 KOL # 
A FROGR OAM FERFORMS KOLMOGOROV —-SMHTRMOV 
A SOORPVMESS—-OF-FIT TEST FOR LOGHORMKMAL 
A ASSUMPTION , K MUST BE GATA ARRAT 
q ORDER DATA AND FIRE SAMPLE MEA ANE VARTAHCE OF LOG TATA 
KEXP A] 
MOM@K 
A FUT TATA TH STANHIARD rHOrRMAL CASE 
Z¢((@K)-KE)=+5 
he POS 
PFE e(.)/O 
IeIIeQ 
eS ,;irtelTi+i 
aL 4xy (IIo) 
q FIHD F VALUE FOR EVERY DATA FOINTS AD FUT Ih ARRAT 
LysiIeI+ 
LAX YC LID 
INVVY Felxo 
sie xX it (ea>s 
oe | 
Eo PPE errr, F 
aL 3 
L4rkKe(y(t-1)) =e 
K&eOQy, K 
CCH (PK) -(PFFF ) 
Che(1t,CC)pFi 
eCe( ,Cry) 
PF ErFFF, OC 
A FI¢ND MAXIMUM ARKSOLUTE DEVIATION ANI FRIMNT AS TEST STATISTIC 
KKe | (CF rPRF-—K ) 
LeKKE PKK J 


) 
OOl 
0) 


Ora On0 msyo-ca bolhore OOO NIO-OIS Whore Ow WII Dib WL 


MMT AGM MAIN AMMAN MMA AMM MMMMmeN 
OL OOP pop hop POP Popo hl i ee be Re re hg 00 SID Li Gihar 
LIL pp LL tL Lt be 


Bee or steatisric ="',¢elf jj 
TVEO.8856=¢M20.5) 
'TARULATEDR VALUE =',eT¥ 


Vv 


¢ 


KOL SET} 
TEST STATISTIC =9,444 


TABULATED VALUE =0,135 


+ Py 


e 102 





Cee REE BY 
VT W KRsALsA2 


3 AS; 
FRODUCES W TEST 


ee ATS TUL 
Tian 


S23 
STATISTIC FOR Ww 


A 

A 

q X= VECTOR THAT IS TESTING 

aA % AND Y MUST BE SAME DIMEHSTONS 


q ORDER DATA 

reek J 

q CALCULATE S22 TERM 
ALe+/ax 


moe (iawn) 

SPeAD-AZ 
q CALCULATE & TERM 

K @ | +3 

LX & @ on 

Lb eer pes 7 

Bet /Vx( (KPLL KJ -—-CKALZ) ) 
Ace teh We TEST STATISTIC AHD FRINT 
Ze (Fx2)+S2 

Steve at StemtriSTtTic = ',¢E 


AMIN IMM MAMMAMMmnAMmnaM:eM, 
Nhe eR RRR Ree ONG On > hk 


mM oOo NO ob Poe OWILILILI DP bIbibil 


bed EJ Los, JLo EEE EL Lk 


340) ) 
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T=eVeGCTCrR FOR W TEST LEHGTH OF ARGUME?HT 
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79 GAMMA ([OI¢9 


9 A GAMMA KovsTp tees age 

Cid 9a INC, GAMMA FHC ,-SHAPE=AQ)9, VECTOR ARG, X30, SCALE FARAM, =} 

C2) 9 GIVES THE LEFT TAIL AREA FOR GIVEN A GHD x 

C3] 9g @ IS SHAPE PARAMETER, XK IS DATA VALUE 

C4] ae. A 

ina J a Ri 

C6] 3OXY.L=A 

ae FEC PRE pX) FO 

C3] Bean eH pve fi 

Cae SL XYO=+/¥ 

meg) 6th S/R 

me el 7eOel / 7) SA+1 

ae J Be tnx, 1 2), Ot yl 

Meee 2(it3+GbcC)x11t<a 

CL4ad FPOV/YVPR JC RX CTRADS CLAD ART 

egg = fk 

VP de rey xex(x/( 1 xe—-1SO-{1 )2,2A~1{ Bj) xCTRA-A) SL, AA 
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